Two-Dimensional Modeling of Soft Ferromagnetic Films 
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We examine the response of a soft ferromagnetic film to an in-plane applied magnetic field. Our 
theory, based on asymptotic analysis of the micromagnetic energy in the thin-film limit, proceeds in 
two steps: first we determine the magnetic charge density by solving a convex variational problem; 
then we construct an associated magnetization field using a robust numerical method. Experimental 
results show good agreement with the theory. Our analysis is consistent with prior work by van den 
Berg and by Bryant and Suhl, but it goes much further; in particular it applies even for large fields 
which penetrate the sample. 



Soft ferromagnetic films are of great interest both for 
applications and as a model physical system. Their sen- 
sitive response to applied magnetic fields makes them 
useful for the design of many devices, including sensors 
and magnetoelectronic memory elements Therefore 
soft thin films have been the object of much experimental 
and computational study Q . Their relatively simple do- 
main structures and significant hysteresis make such films 
a convenient paradigm for analyzing the microstructural 
origin of magnetic hysteresis || . 

Most current modeling of soft thin films is based on 
direct micromagnetic simulation [Q. This is demanding 
due to the long-range nature of dipolar interactions, and 
the necessity of resolving several small length scales si- 
multaneously. Numerical simulation is surely the right 
tool for the quantitative study of hysteresis and dynamic 
switching Q . However it is natural to seek a more analyt- 
ical understanding of the equilibrium configurations. The 
origin of domain patterns is intuitively clear: they arise 
through a competition between the magnetostatic effects 
(which favor pole-free in-plane magnetization) and the 
applied field (which tends to align the magnetization) . A 
2D model based on this intuition was developed by van 
den Berg ](| in the absence of an applied field, and ex- 
tended by Bryant and Suhl Q| to the case of a sufficiently 
weak in-plane applied field, tn van den Berg's model 
(VDBM) magnetic domain patterns are represented us- 
ing 2D, unit-length, divergence-free vector fields, deter- 
mined using the method of characteristics; the caustics 
where characteristics meet are domain walls. In Bryant 
and Suhl's model (BSM) , the presence of a weak applied 
field is accounted for through an electrostatic analogy: 
the "charges" associated with the magnetic domain pat- 
tern should be such as to expel the applied field from 
the interior of the sample, as occurs in an electrical con- 
ductor. The domain patterns predicted by BSM have 
been observed experimentally || . The electrostatic anal- 
ogy is restricted, however, to sufficiently small applied 



fields: since the magnetization vector has a constrained 
magnitude, the field generated by its divergence cannot 
be arbitrarily large. Therefore the BSM breaks down at 
a critical field strength beyond which the external field 
penetrates the sample. 

This Letter extends and clarifies the models of van den 
Berg and of Bryant and Suhl. Our extension is two-fold: 
we permit large applied fields which penetrate the sam- 
ple, and we replace the method of characteristics with a 
robust numerical scheme. Our clarification is also two- 
fold: we identify the regime in which these 2D models 
are valid, and explain their relation to classical micro- 
magnetics. To assess the extended model, we compare 
its predictions to experiments on Permalloy thin film 
elements with square cross-section. The agreement be- 
tween theory and experiment is remarkable, even in the 
field penetration regime. At the heart of our approach 
is an asymptotic analysis of the micromagnetic energy 
in the thin-film limit. The lowest-order terms lead to 
constraints such as = 0, while the second-order term 
sets the charge density. Wall energies and anisotropy con- 
tribute only at higher order. The higher-order terms are 
not irrelevant: they are the source of magnetic hysteresis. 
Our analysis indicates, however, that certain quantities 
should have little or no hysteresis - namely the charge 
density, the region of field penetration, and the magneti- 
zation in the penetrated region. 

The free-energy functional of micromagnetics in units 
of J 2 £ 3 /2^ is 

E d {m) = (nd) 2 / \\7m\ 2 dx + Q ip(m)dx 
Ja d Jn d 

+ I \h d \ 2 dx-2 I h' e -mdx. (I) 
Jr 3 Jn d 

Here m is the magnetization (in units of the saturation 
magnetization J s ), a unit vector field defined on the film 
fid with cross section v and thickness d, where all lengths 



1 



are measured in units of a typical lateral dimension L (the 
diameter for ui a circle, the edge-length for uj a square). 
Moreover, k is the ratio between Bloch line width Dbl 
and the film thickness, where Dbl = {^HqA/ J 2 ) 5 , with 
A the exchange constant, measures the strength of the 
exchange energy relative to that of dipolar interactions; 
Q is the quality factor measuring the relative strength of 
the magnetic anisotropy (p; hd is the stray field in units 
of J s //io, whose norm squared gives the magnetostatic 
energy density; h' e is the applied field in units of J s / Ho, 
which we assume to be uniform and parallel to the film's 
cross section. In what follows, a prime will always denote 
a two-dimensional field or operator. 

For d <C 1 a hierarchical structure emerges in the en- 
ergy landscape of (|l|), see Table |. Variations of m of 
order 1 along the thickness direction X3 give rise to an 
exchange energy per unit area (of the cross section) of 
order n 2 d. An out-of-plane component 1TI3 of order one 
determines a magnetostatic contribution per unit area of 
order d. The component of the in-plane magnetization 
in' orthogonal to the lateral boundary dui of the film's 
cross section u leads to a magnetostatic contribution of 
order d 2 In ^ per unit length. The same mechanism pe- 
nalizes jumps [m' ■ v'\ of the normal component of the 
magnetization across a line of discontinuity of ml with 
normal v' . These lines of discontinuity arise by approxi- 
mating domain walls as sharp interfaces. At order d 2 we 
find the magnetostatic energy per unit area due to sur- 
face "charges" proportional to the in-plane divergence 
div'm'. Finally, the energy per unit length of a Neel or 
asymmetric Bloch wall and the energy of a single vortex 
are indicated in the table. In the regime 
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the highest-order terms penalizing 7713, and [m' • v'\ 
become hard constraints, while the energetic cost of 
anisotropy, of the wall type of minimal energy 9, and 
of vortices become higher-order terms. The energy is 
thus determined, at principal order, by the competition 
between the aligning effect of H' e and the demagnetizing 
effects due to div'm'. 

In view of this separation of energy scales in the regime 
(||) , we propose the following reduced theory. We call an 
in-plane vector field m'(x') on uj "regular" if it satisfies 
[m! ■ v'\ — across all possible discontinuity lines and at 
dui . Our reduced theory states that the magnetization 
m'{x') minimizes 



E(m') = / \H d \ 2 dx-2 / H' e -m'dx', 

where Hd(x) = — \7U is determined by 

V 2 U = in R 3 outside of uj, 
dU' 
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among all regular in-plane vector fields ml of unit length 

|m'| = I in uj. (4) 

Our formula for the induced field Hd is naturally con- 
sistent with that commonly used for 2D micromagnetic 
simulations Eoj. 

We now make two crucial observations. The first is 
that the functional E depends on m! only via the surface 
charge a = —div'm', and it is strictly convex in a. In- 
deed, f R3 \Hd\ 2 dx is a quadratic functional of a and an 
integration by parts shows that J H' e ■ ml dx' is a linear 
functional of a. The second observation is that the set 
of regular in-plane vector fields of unit length and with 
given surface charge is large in the following sense: For 
any regular m' of at most unit length, that is 



< 1 in uj, 



(5) 



there exist many regular ml of unit length with the same 
surface charge: div m' = div m' . Indeed, we may write 
m' = V ^ + m' where V^V' = (—dip/dx2,dip/dxi) and 
the continuous function ijj[x') on uj solves the boundary 
value problem 



if) = on duj. 



(6) 
(7) 



m on uj , 



Condition (g) ensures the solvability of this boundary 
value problem. One can generate many solutions by im- 
posing the additional condition ip = on an arbitrary 
curve contained in uj. 

These observations have two important consequences. 
First, the minimizer of the reduced energy E is not 
uniquely determined. Indeed, according to our first ob- 
servation, E depends only on the surface charge, and 
according to our second observation, a regular in-plane 
vector field of unit length is not uniquely determined by 
its surface charge. 

The second consequence is that the surface charge and 
thus the stray field are uniquely determined. Indeed, ac- 
cording to our first observation, £ is a strictly convex 
function of the surface charge, and according to our sec- 
ond observation, the set of surface charges which can be 
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generated by regular in-plane vector fields of unit length 
is convex. (This is true despite the fact that the set of 
regular in-plane vector fields with unit length is not con- 
vex.) 

Any minimizer ml of (||,(i]) satisfies the Euler-Lagrange 
equation 



H = 0.00 



H = 0.35 



// = 0.57 



H' = Ai 



(8) 



where A(x') is the Lagrange multiplier associated with 
the pointwise constraint (|4|). Since Hd is uniquely deter- 
mined, the region {H' d + H' e ^ 0} of u> where the external 
field is not expelled from the sample is uniquely deter- 
mined. Within this penetrated region, ml is also uniquely 
determined in view of (||). 

There is a finite critical field strength H cr it , in the fol- 
lowing sense: when the applied field is subcritical A = 
and the field is completely expelled from the sample, 
whereas when it is supercritical A is nonzero somewhere 
and the field penetrates in that part of the sample. The 
critical field strength depends on the geometry of uj — 
for a circular disk of diameter 1, its value is 1. Further 
analysis indicates that there can be no walls (discontinu- 
ity lines of ml) in the penetrated region. Moreover the 
penetrated region must meet the boundary of u>. 

To derive quantitative predictions from our reduced 
model (M^), we proceed in two steps. The first step min- 
imizes (0) among all regular in-plane vector fields m' of 
length less than or equal to 1. Recall that replacing (^) 
by ^ does not change the minimum energy; therefore 
the m' obtained this way has the correct reduced energy, 
though it typically violates (Q) . The second step postpro- 
cesses m' by solving (^,^) to obtain another minimizer 
ml of unit length. This ml is the desired magnetization. 

The first step is a convex (though degenerate) vari- 
ational problem. We solve it using an interior point 
method : the convex constraint is enforced by adding 
to the physical energy E a small multiple t of a self- 
concordant barrier B. The unique stationary point of the 
strictly convex E + tB is computed by Newton's method; 
it serves as an initial guess for the minimizer of E + t'B, 
where t' < t. The parameter t is slowly decreased by 
multiplicative increments. Within Newton's method, the 
Hessian of E + tB is inverted by a preconditioned con- 
jugate gradient method. The magnetostatic part of the 
Hessian is evaluated with the help of FFT. This is a ro- 
bust procedure. 

For the second step, we recall that the solution of 
(^,0) is not unique. However there is a special solution 
ip, known as the "viscosity solution", which has special 
mathematical properties JT2| ]. It is robust and can be 
computed efficiently using the "level set method" [ p"3| . 
This is what we compute. 

Our numerical scheme selects — automatically and ro- 
bustly — one of the many minimizers ml . The selection 
principle implicit in this scheme is the same as the one 




FIG. 1. Predictions of the theory: gray-scale plots of the 
vertical component of magnetization. 

proposed by Bryant and Suhl. It appears to pick a min- 
imizer with as few walls as possible. Thus it is not un- 
like the more physical selection mechanism of minimizing 
wall energy, represented as a higher-order correction to 
(I) & 

Figure tt\ shows the predictions of our numerical scheme 
for a square film of edge-length one, subject to a mono- 
tonically increasing field applied along the diagonal. To 
check our predictions, we have observed the response of 
two ac-demagnetized Permalloy (Ni 81 Fei 9 , J s — 1.0 T) 
square samples of edge lengths L = 30 and 60 /im and 
thicknesses D = 40 and 230 nm, respectively, in a dig- 
itally enhanced Kerr microscope. The observed domain 
patterns are given in Figures ||, |^ where the field intensity 
h e , measured in Tesla, is scaled according to 



L hp 
H = -. 

DJ S 



(9) 



Figure [| examines more closely the predictions of our 
theory for \H' e \ close to H cr i t . We have superimposed on 
each gray-scale plot the level curves of the potential v of 
the penetrated field, defined by 

-V« = H' d + H' e . 

Regions where the field lines concentrate are regions 
where Vu ^ 0, i.e., where the external field has pene- 
trated the sample. Within them, (^) implies that ml is 
parallel to Vv . Our theory predicts that the penetrated 
region must meet the boundary of the sample and that ml 
can have no walls in the penetrated region. The pictures 
confirm this, and show quite clearly that two apparently 
independent phenomena - the expulsion of the domain 
walls from the interior of the sample and the penetration 
of the external field - are in fact two manifestations of 
the same event. 

In summary, our model describes the response of a soft 
ferromagnetic thin film to an applied magnetic field. It 
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FIG. 2. Permalloy films: L = 60 fj,m, D = 230 nm. 

determines the micromagnetic energy to principal order, 
and certain associated physical quantities that should 
have little or no hysteresis — the charge density, the 
region of field penetration, and the magnetization in the 
penetrated region. In addition our approach provides a 
specific magnetization pattern which is consistent with 
experimental observation and may well be the ground 
state. Of course, the magnetization of a soft thin film is 
not uniquely determined by the applied field: the multi- 
plicity of metastable states is a primary source of hystere- 
sis. Our approach does not provide a model for hysteresis 
or a classification of stable structures — this would seem 
to require analysis of higher-order terms in the micro- 
magnetic energy. 

We thank S. Conti for valuable discussions. 
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